Parisi States in a Heisenberg Spin-Glass Model in Three Dimensions 
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We have studied low-lying metastable states of the ± J Heisenberg model in two (d = 2) and three 
(d — 3) dimensions having developed a hybrid genetic algorithm. We have found a strong evidence 
of the occurrence of the Parisi states in d = 3 but not in d — 2. That is, in L d lattices, there exist 
metastable states with a finite excitation energy of AE ~ O(J) for L —* oo, and energy barriers 
AW between the ground state and those metastable states are AW ~ 0(JL ) with 9 > in d = 3 
but with 6 < in d — 2. We have also found droplet-like excitations, suggesting a mixed scenario of 
the replica-symmetry-breaking picture and the droplet picture recently speculated in the Ising SG 
model. 
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Recently, the Heisenberg spin-glass (HSG) model in 
three dimensions (d = 3) has been attacked a great in- 
terest, because evidence of the occurrence of SG phase 
transition at a finite, no-zero temperature (To 7^ 0) has 
been given in numerical studies in contrary to a com- 
mon belief that any phase transition does not occur with- 
out some anisotropy [1,2]. Kawamura and his coworkers 
took notice of chiralities of the spins and showed that a 
chiral glass (CG) phase transition occurs at Tcg 7^ 0, 
but the spin glass phase is still absent [3-5]. On the 
other hand, Matsubara et al. examined the stiffness at 
T = and T ^ of the ±J HSG model on the L 3 
lattice with open boundaries and suggested that the SG 
phase transition will occur at T$q ~ 0.19 J [6,7]. They 
also obtained almost the same transition temperature by 
using different numerical methods, i.e., an aging effect 
[8] and the divergence of the SG susceptibility [9]. Re- 
cently, Nakamura and Endoh showed that, using non- 
equilibrium relaxation method, the CG phase transition 
and the SG phase transition occur at the same temper- 
ature of T SG = T CG ~ 0.21 J [10]. Quite recently, Lee 
and Young presented the same conclusion using a finite 
size analysis of the correlation length of the spins and 
chiralities [11]. 

An interesting question is, then, the nature of the SG 
phase of the HSG model. In the Ising SG (ISG) models 
in d = 3, two pictures have been extensively discussed: 
the replica-symmetry-breaking (RSB) picture of Parisi 
[12] and the droplet picture of Fisher and Huse [13]. An 
important difference between the RSB and the droplet 
pictures concerns the nature of large scale excitations. 
In the RSB picture, there are many metastable states, 
which involve turning over a finite fraction of the spins 
and which cost only a finite energy even in the thermody- 
namic limit. The energy barriers between each of these 
metastable states and the ground state and also between 
those metastable states are infinite in the thermodynamic 
limit. Hereafter we call those metastable states Parisi 



states. By contrast, in the droplet picture, the lowest ex- 
citation which has linear spatial extent I typically costs 
an energy of order Jl e with 9 > 0. Hence, in the ther- 
modynamic limit, excitations which flip a finite fraction 
of the spins cost an infinite energy. 

In this Letter, having developed a hybrid genetic algo- 
rithm (HGA) for systems with the XY and Heisenberg 
spins, we have studied the ground state and low-lying 
metastable states of the ± J HSG on finite lattices of L d 
(d = 2 and 3). We have found the Parisi states in d = 3 
but not in d = 2. This finding is very important, be- 
cause it gives a strong support of the presence of the SG 
phase in d — 3. We have also found metastable states 
with droplet-like excitations near the Parisi states. We 
hope our findings will help to understand properties of 
the low temperature phase of the HSG model. 

We start with the ±J Heisenberg model on the d di- 
mensional lattice of L d {= N) with periodic boundary 
conditions. The Hamiltonian is described by 



H 



JijSiSj 



(1) 



where Si is the Heisenberg spin of \Si\ = 1, and (ij) runs 
over all nearest-neighbor pairs. The exchange interaction 
Jij takes on either +J or — J with the same probability 
of 1/2. 

Before analyzing the model, we briefly mention prop- 
erties of the Heisenberg SG model. First we should note 
that the model has the SO(3)(= SU(2)/C 2 ) symmetry. 
So the ground state has a two-hold degeneracy. We call 
these two spin configurations as G\ and G2 and describe 
them {Sf 1 } and {Sf 2 }. Note that, if G x is obtained, 
one can readily obtain G2 by reversing all the spins of 
G\, i.e., {Sf 2 } — { — Sf 1 }, and vice versa. The energy 
of the model is unchanged under any uniform rotation of 
the system. Then the spin configuration of the Heisen- 
berg SG model is characterized by a set of pair-spin cor- 
relations {SiSj}. That is, the SG order parameter space 



1 



is constructed by this set and the distance between two 
spin configurations A with {Sf } and B with {Sf } is 
described as 



N 2 



(2) 



Hereafter we take the ground state as the datum point 
and consider the distance AS of the spin configuration 
A from the ground state spin configuration G\ or G2, 
i.e., AS = S(A,Gi) = S{A,G 2 ). Note that, since the 
ground states G\ and Gi share the same point in the SG 
order parameter space, hereafter we omit the subscript. 
Note also that, if {Sf } is independent of {Sf }, AS ~ 

^^{{Sfsff + (Sfsf) 2 ) ^£/co S (0)2dfi = 

■\/2/3 for i — ► 00. 

Usually, the ground state of the HSG is searched by 
using a spin quench (SQ) method [3,14]. However, the 
number Ni of initial spin configurations which are needed 
to get the ground state increases rapidly as the size of the 
lattice increases. Here we develop a hybrid genetic algo- 
rithm (HGA) to search the ground state and low-lying 
metastable states. Starting with a population N p of ran- 
dom spin configurations (parents) {-B;}, new configura- 
tions (offsprings) are generated by recombination of dif- 
ferent parents Bi and B m (I ^ m) (quadruple crossover). 
Then some fraction r of the spins are refreshed (muta- 
tion) . This algorithm is hybridized with a local optimiza- 
tion (the SQ method) of the offsprings. The population 
N p is updated by replacing 3A p /4 parents with higher 
energy by ZN p /A offsprings which are selected in order of 
the lower energy. This procedure is repeated for many 
times (generation number N g ). This algorithm is anal- 
ogous to that used in the Ising SG model [15]. In the 
Heisenberg SG, however, we should pay a special atten- 
tion to optimize the interface energy between two parents 
Bi and B m when one generates the offsprings. This prob- 
lem can be resolved by applying a uniform rotation to all 
the spins of one of the parents. Having used this algo- 
rithm, we have been able to obtain the ground state of 
the HSG model on larger lattices ofl2xl2xl2. The 
HGA also works well in the XY SG models. Details of 
the HGA will be reported elsewhere. 

Once the ground state G with the energy Eq are 
determined, one can obtain metastable states A with 
Ea- We use the distance AS and the excitation en- 
ergy AE(= Ea — Eg) to distinguish different metastable 
states. We may use two different methods for searching 
those states: (i) the usual SQ method and (ii) the HGA. 
Both methods have their own merit. In the former, we 
can get different metastable states. However, we have 
to start with many different initial spin configurations to 
get low- lying metastable states, e.g., N ~ 10 6 for L = 10 
in d — 3. On the other hand, the latter is appropriate 
to search the lowest metastable state in a given range of 



AS. We are enough to give a smaller number of parents, 
e.g., N p ~ 100 for L = 10 in d = 3. In Figs. 1 (a) and 
(b), we present results of distribution of the metastable 
states of typical bond realizations (samples) in d = 3. In 
fact, by using the two methods we can get the same lowest 
metastable state in a given range of AS (here AS > 0.4). 
Usual excitation energies AE's increase with AS* and also 
with L. A remarkable point is that AE's of the low-lying 
metastable states are much smaller than those usual ex- 
citation energies. In particular, we often see metastable 
states with finite AS which have a very small AE (see 
Fig. 1(b)). This result implies the presence of the Parisi 
states. We will discuss this problem considering an en- 
ergy barrier between the ground state and those low-lying 
metastable states on the basis of a domain wall picture. 
Another remarkable point is that, near the ground state, 
there are many metastable states the energy of which in- 
creases as AS increases. This result suggests the occur- 
rence of some localized excitation the excitation energy 
of which increases as the scale of the excitation is in- 
creased [16]. That is, the excitation will be droplet-like. 
The same is true for the excitations in metastable states 
near the Parisi states. We will get back to this point to 
consider properties of the low temperature phase. 
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for AS > 0.4 indicated by the symbol O- It is interesting 
to see that AW's are much higher than AE. 



10' 



J 1 

i (b) L 


= 10 


1 1 1 

3d , : ^mm 


L- 


^ x 






X. 


















6 




V G , 




i i i 


i i 



"I I I I I I I I ] 

0.2 0.4 0.6 0.8 

AS 



FIG. 1. Distributions of metastable states (dots) in the 
(AS, AE) plane for typical samples with (a) L — 6 and (b) 
L — 10 in d = 3 obtained by using the SQ method with 
N ~ 2 x 10 4 (L = 6) and 2 x 10 6 (L = 10). O indicates the 
lowest metastable state for AS* > 0.4 obtained by using the 
HGA and x denote the domain wall states described in the 
text. G indicates the position of the ground state. 

Suppose that, in addition to the ground state G with 
Eg, some low- lying metastable state A is given. The 
domain wall energy AW between G and A may be es- 
timated by the following procedure, (i) We divide the 
L d lattice into two parts Hi and H2, which are com- 
posed of \(L + l)/2] and [L/2] layers, respectively, where 
[x] means the largest integer which does not exceed x. 
(ii) Fill Hi with the spins of G and H2 with the spins 
of A, and apply a uniform rotation to all the spins on 
H2 to minimize the interface energy, (iii) For each of 
H\ and H2, fix all the spins on the middle layer. Un- 
der this restriction, the SQ method is applied to get the 
spin configuration W which gives the minimum energy 
E\y of the L d lattice, (iv) We consider W a domain 
wall state when S(A,G)/4 < S(W,G) < S(A,G) and 
S(A, G)/4 < S(W, A) < S(A, G) are satisfied, and define 
the domain wall energy AW = Ew — Eq. The num- 
ber of possible divisions of the L d lattice into Hi and 
H2 is d x L and the chirality freedom of 2 exists. So 
we repeat this procedure 2dL times. In Figs. 1(a) and 
(b), we have added AW's of these domain wall states be- 
tween the ground state and the lowest metastable state 



Now we consider the domain wall energies between 
the ground state and the lowest metastable state in a 
given range of AS for different samples. Here we con- 
sider the case 0.4 < A 5 < 0.6, because we arc in- 
terested in large scale excitations. We denote the ex- 
citation energy of the lowest metastable state in this 
range as AEq. Our attention is paid to the minimum 
domain wall energy AW m i n (= min( AW's )) for each 
of those samples. The calculation has been performed 
in d — 2 and d = 3 by using the HGA. The linear 
sizes of the lattice are L = 10^24in<i = 2 and 
L = 6^11inrf = 3, and the numbers of the sam- 
ples are N s = 1024 in both d = 2 and d = 3 except 
for the largest lattices (iV s = 512 for L = 24 (d = 2), 
and N s = 256 for L = 11 (d = 3) ). The following 
numbers of N p and N g are chosen with a common muta- 
tion ratio r = 0.4. In d = 2, N p = 16, 32, 64, 128, 256 for 
L = 10, 12, 16, 20, 24, respectively, and N g = 5 for L < 16 
and N g = 16, 32 for L = 20, 24, respectively. In d = 3, 
N p = 16, 32, 64, 128, 256, 512 for L = 6,7, 8, 9, 10, 11, re- 
spectively, and N g — 5 for L < 8 and N g = 8, 16, 32 for 
L = 9, 10, 11, respectively. Note that, metastable states 
are not always found in this range of AS by using the 
HGA with r = 0.4 [17], especally in smaller lattices. In 
Fig. 2, we plot (Ai?o>AW m i n ) for different samples in 
d = 3. A remarkable point is that, as L is increased, 
the distribution of AW m i n considerably spreads to the 
higher energy side, while the distribution of AE slightly 
shifts to the lower energy side. Note that, in d = 2, the 
distributions in both AEq and AWmin shift to the lower 
energy side with increasing L. 
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FIG. 2. The lowest excitation energy AEo vs the minimum 
domain wall energy AWmin of different samples for (a) L — 6 
and (b) L — 10 in d — 3. Here N s is the number of the sam- 
ples in which some metastable states with 0.4 < AS < 0.6 
are obtained among 1024 samples. 



We calculate average values (AEo) and (AW m i n ) over 
different samples and show them in Figs. 3 and 4 for 
d = 2 and d = 3, respectively, as functions of L. In 
d = 2, in fact, (AEo) and (AW m i n ) decrease with in- 
creasing L. These results clearly reveal the absence of 
the SG phase. By contrast, in d = 3, (AEo) decreases 
slightly, and (AW m i n ) increases with L. We could fit data 
for larger L as (AW min ) oc JL 9 with 9 = 0.53 ± 0.08. 
These results strongly suggest the presence of the Parisi 
states. That is, in the thermodynamic limit, there are 
metastable states which have a finite excitation energy 
AEo, probably AE n <C J, and which are separated from 
the ground state with an infinite energy barrier. The ex- 
ponent 9 is the measure of the domain wall height. It is 
interesting to find that the value of 9 ~ 0.53 is compati- 
ble with the stiffness exponent of 9 = 0.4 <~ 0.8 estimated 
recently [6,7]. 



FIG. 4. Average values of the lowest excitation energy 
(AEo) and the lowest domain wall energy (AW m i n ) in d = 3 
as functions of the linear lattice size L. 

In summary, we have studied low-lying metastable 
states of the ±J Heisenbcrg model in two (d = 2) and 
three (d = 3) dimensions having developed a hybrid ge- 
netic algorithm. We have found the Parisi states in d = 3. 
Since the Parisi states occur in the SG phase, we suggest 
that the SG phase really occurs in the HSG in d = 3. We 
have also found metastable states with droplet-like exci- 
tations near the Parisi states. So the SG phase of this 
HSG model will be characterized by a mixed scenario of 
the RSB picture and the droplet picture analogous to 
one [18] that was speculated previously in the Ising SG 
model [19]. In this scenario, the system at T — is lo- 
cated in one of the Parisi states and dynamical properties 
of the system at low temperatures may be described by 
using the droplet-like excitations thermally activated in 
that Parisi state. In the Ising model, aging dynamics 
of the SG model in d = 3 has been reported to be well 
described by using the droplet picture [20], although the 
model has large scale excitations with finite excitation 
energy [21,22]. We hope that the present findings stimu- 
late studies of aging dynamics of the HSG as well as the 
equilibrium properties of the real SG systems. 
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FIG. 3. Average values of the lowest excitation energy 
(AEo) and the lowest domain wall energy (AW m i n ) in d = 2 
as functions of the linear lattice size L. 
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